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Abstract 

We discuss the generation of large-scale magnetic fields due to the breaking of the conformal 
invariance in the electromagnetic field through the CPT-even dimension-six Chern-Simons-like 
effective interaction with a fermion current in inflationary cosmology. It is shown that the magnetic 
fields on IMpc scale with the field strength of ~ 10 -9 G at the present time can be generated even 
for the scale of the effective interaction being the Planck scale. 

PACS numbers: 98.80. Cq, 98.62.En 
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I. INTRODUCTION 



It is observationally known that there exist magnetic fields with the field strength 10~ 7 - 
10 _6 G on lOkpc-lMpc scale in clusters of galaxies as well as ~ 10 _6 G on 1 — lOkpc scale 
in galaxies of all types and at cosmological distances (for reviews, see [l]). However, the 
origin of the cosmic magnetic fields, in particular the large-scale magnetic fields in clusters 
of galaxies, is not well understood yet. Although the dynamo amplification mechanism 
amplifies very weak seed magnetic fields up to ~ 10 _6 G, its effectiveness in galaxies at 
high redshifts and clusters of galaxies is still unsatisfactory. Furthermore, it is difficult 
for astrophysical processes 0, 4], cosmological phase transitions [5] and primordial density 
perturbations before or at the epoch of recombination [6] to generate the magnetic fields on 
megaparsec scales with sufficient field strengths to account for the magnetic fields observed 
in galaxies and clusters of galaxies without dynamo amplification mechanism. 

The most natural origin of the large-scale magnetic fields is from electromagnetic quan- 
tum fluctuations existed at the inflationary stage 3). This is because inflation naturally 
produces effects on very large scales, larger than the Hubble horizon, starting from micro- 
physical processes operating on a causally connected volume. In the Friedmann-Robertson- 
Walker (FRW) spacetime, the metric is conformally flat, while the ordinary Maxwell theory 
is conformally invariant. It is clear that the conformal invariance must have been broken 
at the inflationary stage in order that electromagnetic quantum fluctuations can be in- 
duced during inflation [8j. Note that this does not apply when the FRW background has 
nonzero spatial curvature [9]. Various breaking mechanisms of the conformal invariance in 
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22j | . and cosmic defects |23j, have been proposed. 



Recently, an C PT-even dimension-six Chern-Simons-like effective interaction between a 
fermion current and the electromagnetic field in inflationary cosmology was proposed to in- 



duce the cosmological birefringence 24], |25( as well as baryon number asymmetry [26| . The 



electromagnetic quantum fluctuations can be generated during inflation as the interaction 
breaks not only the Lorentz invariance but also the conformal invariance of the electromag- 
netic field. In the present paper, we study the generation of large-scale magnetic fields due 
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to this effective interaction. 

This paper is organized as follows. In Sec. II we describe our model and derive equations 
of motion for the U(l) gauge field. In Sec. Ill we consider the evolution of the U(l) gauge 
field and estimate the present strength of the large-scale magnetic fields. Finally, Sec. IV is 
devoted to a conclusion. We use units of k-Q = c = H = 1 and adopt Heaviside-Lorentz units 
of electromagnetism. 



II. MODEL 



The action for the Maxwell theory with the C PT-even dimension-six Chern-Simons- 
effective interaction between a fermion current and the electromagnetic field (A^ 
is given by 
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where g is the determinant of the metric tensor g^, F^ u = d^A u —d u A^ is the electromagnetic 
field strength tensor, F^ = [1/ (2- s /=g))e tMUpa F pa is the dual of F^ with e^ vpa being the 
Levi-Civita tensor normalized by e 0123 = +1, /3 is a dimensionless coupling parameter, 
and M = A /Air with A being the scale of the effective interaction. The Chern-Simons- 
like effective interaction might originate from the low energy effective theory in superstring 
theory [24]. 

From the action in Eq. (12.11) . the equation of motion for the electromagnetic field can be 
derived as 
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We take the flat FRW space-time with the metric, ds 2 = g^dx^dx" = —dt 2 +a 2 (t)dx 2 , where 



a(t) is t 



form 



re scale factor. The fermion current to a comoving observer has the following 
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where and are the number densities of the fermion ip and antifermion if>, respectively. 
In the FRW background, for the Coulomb gauge of A (t, x) = and the case of djA^it, x) = 



0, Eq. (|22D becomes 

Ai(t, x) + HAi(t, x) - \djdjAiit, x) + 2j^na- 1 e ijk d j A k (t, x) = , (2.4) 

where a dot denotes a time derivative, H = a/a is the Hubble parameter, and e^-fe is the 
totally antisymmetric tensor (ei23 = +1). 



III. LARGE-SCALE MAGNETIC FIELDS 



A. Evolution of the U(l) gauge field 



We now consider the case in which a slow-roll exponential inflation is realized with the 
scale factor a(t) given by a(t) = a\ exp [H- m f(t — t\)] , where a\ is the scale factor at the time 
t\ when a comoving wavelength 2ir/k of the U(l) gauge field first crosses outside the horizon 
during inflation, k/(aiH in f) = 1, and H in f is the Hubble constant at the inflationary stage. 
From the quantization of the U(l) gauge field A^t, x), we obtain the expression for Ai(t, x) 
as 

d 3 k 
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b(k)Ai(t, k)e ikx + b\k)A*{t, fc)e~ ifc ' a 



(3.1) 



where k is the comoving wave number, k denotes its amplitude |fe|, and b(k) and tf(k) are the 
annihilation and creation operators which satisfy b(k), b^(k') = § 3 (k — k') and others = 0. 
In what follows, we choose the x 3 axis to lie along the spatial momentum direction k and 
denote the transverse directions x 1 with I = 1,2. We use circular polarizations expressed 
by the combination of linear polarizations as A±{k, t) = Ai(k, t) ±iA 2 (k, t). From Eq. (12.41) . 
we find that 
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Afc(M) + H ini A±(k,t) + 
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A±(k,t) = 0. 



(3.2) 



Unfortunately, it is impossible to obtain the analytic solution of Eq. (13.21) for a generic 
evolution of the fermion number density n at the inflationary stage. If n evolves as n oc 
a oc —1], where rj — J dt/a(t) is conformal time, or n oc a -1 , analytic solutions for Eq. 
( 13. 2p can be derived. We will investigate these cases in the next subsection. Here we will 
numerically solve Eq. (13.21) at the inflationary stage. We assume that the initial amplitudes 
of A + (k,t) and A_(k,t) are the same with the time t\ as the initial time. During inflation 



(ti < t < tn, where t-& is the end of inflation), the amplitudes of A±(k, t) are expressed 
as A±(k,t) = C±(k,t)A±(k,ti), where C±(k,t) are obtained by numerical calculations with 
C±(k, tx) = 1. We note that C±(k, t) has the ^-dependence. Hence, in numerical calculations 
we chose a comoving scale L = 2^/k. By requiring that the vacuum should reduce to the 
one in Minkowski spacetime in the short-wavelength limit, we have | (fc, ) | = \ j\fTk and 
|At(MOI =H iv t/V2k. 

For n = na~ 3 , the numerical results for the evolutions of C + (k, t) and C_(/c, t) are shown 
in Figs. 1 and 2, respectively, where n is a constant, £ n = (n/ [cur 3 ]) 1/2 / [Mj [GeV]) = 
4.53 x 10" 44 , H ini = 10 10 GeV, (3 = 1.0 and a comoving scale L = 2ir/k = IMpc. By using 
the five year Wilkinson Microwave Anisotropy Probe (WMAP) data on the anisotropy of 
the cosmic microwave background (CMB) radiation [^tJ, one gets that H m f < 6.0 x 10 14 GeV 
from tensor perturbations [28(. As shown later, the large-scale magnetic fields with the 
sufficient amplitude at the present time can be generated. Here, we have used k/a = 
exp [— H m i {t — ti)} H in i and a\ = k/H in f, and calculated Eq. ( 13.21) numerically at t = t\ = 
Hrj with C±(Mi) = 1-0 and C±(Mi) = H in{ . 

The numerical results in Figs. 1 and 2 are understood as follows. From Eq. ( \3.2\\ for 
n = na~ m with m being an integer, the equation of C±(k, t) is given by 

ci (k,t) + c ± (k,i) 

+ exp[-2(t-fi)] {lTJexp[-{m-l)(i-i 1 )]}C±(k,i) =0, (3.3) 

where 

t = H ini t, t\ = H iDi ti, (3.4) 

and the prime denotes the derivative with respect to t. We now consider the case of m = 3 
shown in Figs. 1 and 2. At the early stage of inflation, the second term in the braces { } 
on the left-handed side of Eq. (I3.3P is much larger than the first one. This means that 
the conformal invariance of the electromagnetic field is broken, and hence the amplification 
of C± (k,t) can be realized. However, as t becomes much larger than t 1; the second term 
becomes much smaller than the first one due to the existence of the exponential suppression 
term, exp [—2 (t — . After about ten Hubble expansion times, the second term become 
negligible comparing with the first one, and the above equation becomes almost equal to 
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FIG. 1: Evolution of C+(k,t) for n = na~ 3 , £ n = (n/ [cm- 3 ]) 1/2 / (M/ [GeV]) = 4.53 x 10" 44 , 
^inf = 10 10 GeV, P = 1.0 and a comoving scale L = lixjk = IMpc. 




FIG. 2: Evolution of C_(fc,t) for n 



na 



-3 



Legend is the same as Fig. 1. 



the equation of the ordinary Maxwell theory, which can be rewritten to the following form 
by replacing the independent variable t with the conformal time 77: 

d 2 A±(k,r]) 



drj 2 



+ k 2 A ± (k : rj) = 0. 



(3.5) 
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The solution of Eq. (13.51) is given by A±(k, rj) — 1/ y/2ke~ lkri . Thus, the absolute value of the 
amplitude is constant as \A±(k } 77) | = 1/V2k. Consequently, the solutions of C±(k, t) become 
asymptotically constant and the behaviors of C±(k,t) in Figs. 1 and 2 are reasonable. 



B. Strength of the large-scale magnetic fields 

Next, we estimate the present strength of the large-scale magnetic fields. The proper 
magnetic fields are given by [lo| 

Bf TopeT (t, x) = a~ x Bi{t, x) = a~ 2 e ijk djA k (t, x) , (3.6) 

where Bi(t, x) is the comoving magnetic field. From Eq. ( 13. 61) . the energy density in Fourier 
space is 

PB (k,t) = \ [ii?r pcr (M)i 2 + i^ roper (M)i 2 ] , (3.7) 

\Br pev (k,t)\ 2 = -J-) 2 \A ± (k } t)\\ (3.8) 

where B± oper (k, t) = B propcT (k,t) ± iB propeT (k,t). Multiplying p B (k,t) by the phase-space 

density of 4^k 3 /(2n) 3 , we get the energy density of the proper magnetic field as 

P r i 
p B (L,t) = — [\B pmpe \k,t)\ 2 + \B p ^(k,t)\ 2 \ , (3.9) 

in the position space on L. Using Eqs. (I3.7H ( T379T) . we find 

1 /P 4 



PB(L,t) = — - l(k,t), 
8n z \a J 

I(k,t) = \C + (k,t)\ 2 + \C4k,t)\ 2 , (3.10) 



where X(fc, t) corresponds to the amplification factor at the inflationary stage. Here, we 
concentrate on the case in which after inflation the universe is reheated immediately at 
t = tR. The conductivity of the universe a c is negligibly small during inflation because there 
are few charged particles at that time. After reheating, a number of charged particles are 
produced so that the conductivity immediately jumps to a large value: a c ^> H. For a large 



enough a c , magnetic fields evolve in proportion to a 2 (t) 10]. From B(L,t) = ^2p B (L, t) 
and Eq. (I3.10p . we find that the present strength of the magnetic fields is 

B(L, t ) = 8.2 x 10 18 exp (-2N) (pfy) (^) V 7 ^^) [G] , (3.11) 



where a R /a = (g R /3.91) 1//3 T 70 /T R with T R being the reheating temperature and 



291 ]. a R and ao(— 1) 



70/ 

T 7 o (= 2.73[K]) the present temperature of the CMB radiation 
are the values of a at t = t R and the present time to, and N is the number of e- 
folds between the time t\ and t R , given by N = 45 + ln(L/[Mpc]) + In J, where J = 
[30/(tt 2 ^ r )] 1/12 pr 1/4 / (l0 38 / 3 [GeV]) , g R (^ 100) is the total number of degrees of freedom 
for relativistic particles at the reheating epoch, and p R = (7r 2 /30) g R T R 4 is the energy density 
of radiation at the reheating stage. 

Using Eq. f[3TTj) and if 2 nf = (8tt/3) p R /M Ph where M Pl is the Planck mass, we find that 
for n = na~ 3 , £ n = 4.53 x 10~ 44 , F inf = 10 10 GeV and (3 = 1.0, C + (k,t R ) = 1.9 x 10 48 and 
C_(/c,t R ) = —4.6, and consequently the field strength of the generated magnetic fields on 
IMpc scale at the present time is B (L = lMpc,t ) = 1-1 x 10 _9 G. 

We account for the relation of the amplitude of the magnetic fields and its scale depen- 
dence to the evolution of the fermion number density as n = na~ m . It follows from Eq. (13. 3p 
that when the value of J is large, the breaking magnitude of the conformal invariance of the 
electromagnetic field becomes large and hence the strong magnetic fields can be generated. 
From Eq. ([H]) with a x = k/H inf , we find J = 2/3 (n/M 2 ) (1/H inf ) (k/H in{ )~ m oc ^tf™" 1 *;-™ 
When m > 0, the spectrum of the resultant magnetic fields is red one, and the magnetic 
fields on some given large scale, i.e., small k, become stronger as m is larger. 

For the strength of the primordial magnetic fields, there are constraints from Big Bang 
Nucleosynthesis (BBN) on smaller scales. The limit on the present strength of the magnetic 
fields around the BBN horizon size ~ 9.8 x 10~ 5 /i _1 Mpc is less than 10 _6 G 30j. Here, h 
is the related quantity to the present Hubble parameter as Hq = 2.13/i x 10 _42 GeV 29]. 
Note that throughout this paper, we use h = 0.7 [31[. In the case in Figs. 1 and 2, the 
present strength on the BBN horizon scale is 5.9 x 10 _50 G. Clearly, the constraints from 
BBN are satisfied. On the other hand, there exist constraints from the CMB anisotropy 
measurements on larger scales. The result of ~ 10 _9 G in Figs. 1 and 2 on IMpc scale is 
consistent with the recent observational upper bound derived by using the WMAP 5 years 
data [321 ] . Similar bounds have been studied in Ref. 33(. According to Ref. 34j], the limit 
on the current strength on scales larger than the present horizon is less than 4.8 x 10~ 9 G. 
For i n = 1.87 x 10^ 49 , H in{ = 10 10 GeV and p = 1.0, the present strength of the magnetic 
fields on the horizon scale is 2.4 x 10~ 9 G, which is consistent with the above upper limit. In 
this case, the field strength on IMpc is 1.2 x 10~ 57 G. 
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In addition, we study the case in which n evolves as n = na. In this case, we can obtain 
the following analytic solution of Eq. (13. 2p : 



1 , n p -, 1 



where H$p is a z/th-order Hankel function of type 1, and we have taken v > and chosen 
the integral constant so that the vacuum reduces to the one in Minkowski space-time at the 
short- wavelength limit. Being interested in large-scale magnetic fields, we investigate the 
behavior of this solution in the large- wavelength limit. Expanding the Hankel function in 
Eq. (13.131) and taking the first leading order in kj (aH in f), from Eqs. (I3.7p - (I3.9I) we find that 
the present strength of the large-scale magnetic fields is given by 

B ^«- 5 - lxl0l9 S r M(l^j) 2 fe) 2 (^) 5/ " [G| ' (314) 

If v = 2.57 and H m { = 10 14 GeV, the present strength of the large-scale magnetic fields on 
IMpc scale is B(L = lMpc,t ) = 1.4 x 10 -9 G. In this case, £ n = 3.27 x 10 53 . 

It is interesting to note that for n = na -1 , we can also obtain an analytic solution of Eq. 
(13. 2p . In this case, Eq. (13. 2p is rewritten to the following form by replacing the independent 
variable t with the conformal time rj: d 2 A± (k, r/^ / (drj 2 ) + k 2 A± (jc, i]j = 0, where k 2 = 
k 2 =p 2 (P/M 2 ) nk. The solution of this equation is given by A±(k, rf) oc exp (^—ikr^j . This is 
an oscillating solution and its absolute value is constant. Hence A±(k, t) cannot be amplified. 

We remark that if ip is an unknown particle and n = ria ~ 3 . the present number density n 
should be smaller than that of the neutrino 1.1 x 10 2 cm -3 29(. For n/ [cm -3 ] = 3.1 x 10~ 49 , 



1 /2 

which satisfies the above constraint, it follows from £ n = (n/ [cm -3 ]) / [Mj [GeV]) = 
4.53 x 10 -44 that M can be the Planck scale, i.e., M = M P1 = 1.2 x 10 19 GeV. 

In this paper, we treat the fermion ip being relativistic during inflation. The energy 
density of the inflaton is given by p; n f = [3/ (8tt)] H 2 n{ M^ v The energy density of radiation 
at the reheating stage is given by = (7r 2 /30) (?rT r 4 . Since we consider the instantaneous 
reheating in our study, p; n f = pr. The energy density of relativistic fermion is given by 
Pfermion = (7/8) (7r 2 /30) g R T R 4 [29j . Thus, pfermion = (7/8)p inf , namely, the energy density of 
the fermion if> is smaller than that of the inflaton. This is reasonable because in the standard 
inflationary cosmology the potential energy of the inflaton is mainly responsible for inflation. 



Moreover, in this paper we do not study the production mechanism of the fermion ip. We 
assume that it is produced by some other mechanism. 

In addition, we note that the present value of the ratio of the Chern-Simons interaction 
term to the Maxwell one is given by \£cs/£m\ ~ 10/3^, where we have used \d fl A u /A u \ w H . 
Hence, for n = na" 3 , £ n = 4.53 x 10~ 44 and (3 = 1.0, the above ratio is much smaller than 
unity. Thus, the relative contribution of the Chern-Simons interaction term at the present 
time is very small. 



IV. SUMMARY 



In summary, we have studied the generation of the large-scale magnetic fields due to the 
breaking of the conformal invariance in the electromagnetic field through the C PT-even 
dimension-six Chern-Simons-like effective interaction with a fermion current in inflationary 
cosmology. We have found that the magnetic fields on IMpc scale with the present amplitude 
of ~ 10 _9 G can be generated. This strength is enough to explain the magnetic fields observed 
in galaxies and clusters of galaxies through only adiabatic compression without requiring 
any dynamo amplification |7| . If the number density of the fermion ip interacting with the 
electromagnetic field evolves in proportion to a~ 3 (t) during inflation, the scale of the effective 
interaction can be the Planck scale. 



Acknowledgments 

This work is supported in part by the National Science Council of R.O.C. under Grant 
#: NSC-95-2112-M-007-059-MY3 and National Tsing Hua University under Grant #: 
97N2309F1. 



Note added. — After this work was completed, we became aware of a related work 35] 



considering the generation of the primordial magnetic fields during inflation in a Lorentz 



violating theory of Electrodynamics due to a Chern-Simons term coup 



field to an external four- vector, proposed by Carroll, Field and Jackiw 36(. In the scenario 



ing the U(l) gauge 



of Ref . [35[ , during inflation the induced magnetic fields are peaked on much smaller scales 
than the Hubble horizon at that time. To extend the scales to the galactic one at the time 
of protogalactic collapse, the inverse cascade mechanism has to work. On the other hand, 
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in our scenario, during inflation the magnetic fields can be generated on much larger scales 
than the Hubble horizon at that time so that the magnetic fields with sufficient strength 
on the scale of cluster of galaxies at the present time are produced without any secondary 
extension mechanism such as the inverse cascade. This is the advantaged feature of our 
scenario. 
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